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ABSTRACT 

An analysis is made of temperature distributions in a slab of finite 
thickness and infinite extent, one surface of which is perfectly insulated 
and the other surface of which is exposed to a fluid the temperature of 
which varies in a specified way. Two cases are considered. In the first 
case, the temperature of the fluid is suddenly changed and maintained at 
the new value. In the second case the temperature of the fluid increases 
linearly with time. From the solutions of these problems, curve sheets 
are developed which permit evaluating each of eight different quantities 
of engineering significance. These curve sheets cover more cases of 
interest and are easier to use than similar curves which have appeared in 
the report literature. 
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I. INTRODUCTION 



The problem of convective heat transfer by a fluid flowing through a 
thin-wall pipe with a well insulated exterior surface is of considerable 
importance, particularly in connection with assessing the effects of rapid 
operating transients in nuclear power systems. Several years ago contractors 
engaged in this sort of work began using a set of several sheets of curves 
which facilitate these calculations. The writer has been unable to learn 
the original source of these curves or even to find them anywhere in the 
open literature. However, they are to be found in such documents as Refer- 
ences 2 and 3. 

The theory upon which these curves are based is well known and will be 
developed, for reference purposes, in Section II of this thesis. It re- 
places the thin-wall tube by a flat plate which extends infinitely in the 
y and z directions and is of finite thickness a in the x direction. At 
the surface x = 0 it is presumed that there is a perfect insulating layer 
and at the surface x = a the plate is exposed to a fluid whose temperature 
varies with time in a specified manner; the coefficient of surface heat 
transfer between the fluid and the plate is a specified constant h. The 
initial temperature distribution in the plate is given as a function of x. 

The problem is to determine the subsequent temperature distribution in the 
plate as a function of x and the time t. It is further presumed that the 
appropriate physical properties of the plate are constants. 

This problem can be reduced to that of solving a certain linear partial 
differential equation subject to certain initial and boundary conditions. 

To simplify and systematize the analysis it is convenient to introduce and 
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deal with several dimensionless ratios. There is more than one form in 
which the solution may be represented but we will obtain and develop a 
particular and convenient form. 

In the case of a pipe the wall thickness of which is an appreciable 
fraction of its radius, cylindrical geometry should be employed and an addi- 
tional geometric parameter, ratio of wall thickness to radius, must be intro- 
duced. The solution in this case is more complicated than that which will 
be developed here and which forms the basis for the useful curves referred 
to above. However, in most practical cases the wall thickness is indeed 
sufficiently small compared to radius to justify using the flat-plate geo- 
metry which will be employed here. No attempt is made in this thesis to 
explore the differences between results using flat-plate geometry and re- 
sults using cylindrical geometry. 

Two particular cases arise frequently in practice. In each case it is 
assumed that the fluid and the material of the pipe are at the same tempera- 
ture prior to the time t = 0. In one case, referred to as "step", it is 
assumed that the fluid suddenly assumes and maintains a new and different 
temperature starting at time t = 0. In the other case, referred to as 
"ramp* 1 or "linear 11 , it is assumed that starting at time t = 0 the temperature 
of the fluid changes linearly with time. By superposition of results, the 
cases of piecewise linear changes in fluid temperature including sudden 
"jumps" in temperature can be accommodated, if convenient solutions are 
available for the two cases described. 

The useful curves referred to above give significant information 
relating to these two cases. However, those sets of curves which have 
come to hand are considerably less useful than would be desirable because 
they do not provide a sufficiently broad range of significant parameter 
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values to represent physical situations actually encountered and they require 
some rather drastic interpolation and extrapolation. Furthermore, certain 
integral properties of the solutions are useful for certain types of stress 
calculation and the curve sets previously available have not provided these 
properties. It should also be noted that some of the information contained 
In the curves which are developed could also be obtained from the well known 
"Heisler charts,” specifically from Figures 7 and 10 of Reference 4. 

Accordingly, it has been a major purpose of this thesis to reconstruct 
the curve sets, including a greater variation in significant parameters and 
to augment the collection of curves by adding several new curves likely to 
prove useful in engineering practice. The curves developed were checked and 
verified to be correct by comparison to previously published curves. The 
only exceptions to this were the /\T^ and AT^ cases for which there were 
no comparison curves or data available. 

The development, including the differential equation and boundary and 
initial conditions, leading to the production of the desired curves is 
presented in Section II. A. of this thesis. In Section II. B. the theory of 
the step fluid temperature change is developed and that of the linear fluid 
temperature change is developed in Section II. C. A discussion of the compu- 
tational method used to obtain the curve sets is presented in Section II. D. 
Additionally, the computer program has many comments inserted to aid in 
understanding it. A brief discussion of how to use the curve sets is pre- 
sented in the first part of Appendix A, following which the curves them- 
selves are presented. Further instruction is given in References 2 and 3. 

The person who is engaged in heat transfer calculations involving slabs or 
thin wall pipes can easily develop further methods as required. 
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Inasmuch as it is likely that in nuclear power plant calculations there 
may be need to perform calculations relating to convective heat transfer 
of the type described above on a massive scale, involving hundreds of differ 
ent parameter sets and hundreds of transient specifications, it is of in- 
terest to consider ways in which human intervention is not required to "read 
curves# Of course, a computer may be used to calculate each of these possi- 
bly hundreds of thousands of cases using the exact theory developed herein# 
However, since the calculations are, as will be seen, quite demanding, there 
is reason to seek an approximation which will have sufficient accuracy for 
most purposes but which will provide significant speed-up of the calcula- 
tions# For this purpose, in Section III hereof, we seek to determine the 
coefficients of a quickly-evaluated approximation in such a manner that the 
error will be limited to small values over the entire significant ranges 
of all parameters# 

Figure 1 illustrates the problem which will be dealt with# 
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We will find a convenient dimensionless temperature T(x,t)/T^ for the 
two previously mentioned types of fluid temperature change. The step 
temperature change, where = 0 for time t <C 0 and = F = const, for 
t y- 0, will be considered first in Section II. B. The linear temperature 
change, where T f = 0 for t < 0 and T f = Et for t 0 will then be con- 
sidered in Section II. C. Here E = const, is in degrees per unit time. 

Actual solutions will be cast in the dimensionless form T(x,t)/Tj, which will 
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be an explicit function of x or x/a and Fourier number and an implicit 
function of Biot number. 

Each of the two general equations giving T = T(x,t) will then be employed 
to provide eight particular properties most likely to be of use to the 
designer. Curves representing these properties in a convenient way will be 
obtained and exhibited in Appendix A. These results include those which 
have appeared in the curve sets of References 2 and 3 but with a fuller 
coverage of significant parameter values. Additionally some new properties 
are presented that are intended to assist in performing evaluations according 
to the rules of the ANSI B31.7 Code for Nuclear Power Piping and a forth- 
coming edition of Section III, Nuclear Pressure Vessels, of the ASME Boiler 
and Pressure Vessel Code. 
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II. THEORY AND DEVELOPMENT OF EQUATIONS FOR HEAT TRANSFER 



A. GENERAL THEORY 

It is desired to find T(x,t) such that it satisfies the unsteady heat 



conduction equation in one dimension, namely 

„ c^T dT 

dx^ = dx (i) 



The initial and boundary conditions of the problem are: 



1) T(x,0) = 0 (slab initially at uniform temperature throughout, 

and we take this to be zero without loss of 
generality) 



2 ) 

3) 



il i 

ox I x=o 



c>X lx - <X 



0 (perfect insulation at x = 0) 




(continuity of heat transfer at 
fluid-solid interface) 



B. STEP FLUID TEMPERATURE CHANGE 

Here T^ = F = const, for t > 0. The solution to this problem is 
found^ to be 



t( x >* ) - /- g c h £°sCP h 'x)e-" f>h 

F 



( 2 ) 



where : 

a) Ph = 

b) M n Y^/«. 

C) /3h- 



^Schneider, Paul J., Conduction Heat Transfer , p. 250, Addison Wes ley , 

1955. 
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d) C =4 sin M /(2M + sin 2M ) 

n n n n 

e) 2 denotes summation over n = l,2,3,...oo 

f) The numbers are an infinite number of distinct roots as 

seen in Figure 2 below, of the equation 

M tan M = Ng (3) 



where N is the Biot number 
b 

N = ha/k 

-D 




Figure 2 

To verify that the governing differential equation is satisfied we find 
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( 4 ) 
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The governing equation is then 



F£ ChP.'cvsCP.,*)^" 7 '*- F£e n X?c°r (fax) 



e 






But cxMn/o^ Therefore the differential equation is 

satisfied. Next, we consider the first boundary condition which takes 
the form 

/= 2 c h cos(8 h x) 

and we will show that this is true by actually determining the quantities 
C . These are similar to Fourier coefficients which must represent the 
constant value of unity. Even though the argument (f) x ) is not simply 
a multiple of n, the customary technique of evaluation can still be used. 
Multiply through by cos( x)dx and integrate from 0 to a. On the left 



we get 



On the right we get 



X C.vsC&t* " S JJl Hlz 



(ft* 



*2 c h X cos (0 ni x) £ *) J * 



and it is necessary to distinguish between the cases where m = n and m ^ n. 



For m^n, we find . 

. /• „ \ , ^ [ Sth (Mr,-. '/‘I j Sln(Ahy.ifi 1’-,) 

X cosCe*. xj cosCfr*) x l_ - -7 k~7nT l 



/An, CVS _ /lolf/Vj,,, An M r A - Mh M/J 

JL -I 



- «=v. 

which vanishes since M and M are distinct roots of M tan M = N 
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For m = n, we find 
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7. Air, J 



Thus, we have established that 
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57 n Ain = 



o. 

Ain 



r,4 Slh ' 1 Mh l 

n L -2.A'.* J 



so that C =4 sin M /(2M + sin 2M ) 

n n n n 

as previously defined. 

That the second boundary condition is satisfied may be immediately 
verified by inspection of equation (4). Additionally, upon noting that 

#> 57/0 ( fn c\) = 5//) = A <5.5/5 

it is trivial to show that the third boundary condition is satisfied. 

Thus the governing equation and boundary conditions have been shown to be 
satisfied. Equation (2) is interesting and useful ’in its own right. 
However, we now proceed to derive formulas for eight particularly useful 
consequences of the relation. In these formulas the distance x no longer 
appears as a parameter and the dependence upon time is incorporated in the 
value of the Fourier number 



(5) = Fourier number =* 



which is simply related to the exponent in the general formula, viz: 

A= M " 0 

The eight cases we consider and the corresponding formulas are: 



1. Case 1 . The average temperature of the slab is 
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